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ABSTRACT 


Observations  of  winds,  waves,  and  turbulence  at  the  ocean  surface  are  compared  with  several  analytic 
formulations  and  a  numerical  model  for  the  input  of  turbulent  kinetic  energy  by  wave  breaking  and  the 
subsequent  dissipation.  The  observations  are  generally  consistent  with  all  of  the  formulations,  although  some 
differences  are  notable  at  winds  greater  than  15  m  s_1.  The  depth  dependence  of  the  turbulent  dissipation  rate 
beneath  the  waves  is  fit  to  a  decay  scale,  which  is  sensitive  to  the  choice  of  vertical  reference  frame.  In  the 
surface-following  reference  frame,  the  strongest  turbulence  is  isolated  within  a  shallow  region  of  depths  much 
less  than  one  significant  wave  height.  In  a  fixed  reference  frame,  the  strong  turbulence  penetrates  to  depths 
that  are  at  least  half  of  the  significant  wave  height.  This  occurs  because  the  turbulence  of  individual  breakers 
persists  longer  than  the  dominant  period  of  the  waves  and  thus  the  strong  surface  turbulence  is  carried  from 
crest  to  trough  with  the  wave  orbital  motion. 


1.  Introduction 


Gemmrich  and  Farmer  2004;  Gemmrich  2010).  The 
dynamic  balance  assumed  is  that  the  wave  energy  lost 
during  breaking  Sbrk  becomes  a  flux  of  turbulent  ki¬ 
netic  energy  F  into  the  ocean  surface  layer  that  is  dis¬ 
sipated  at  a  rate  e(z),  which  varies  with  depth  z 
beneath  the  surface: 


Wave  breaking  at  the  ocean  surface  limits  wave 
growth  (Melville  1994),  enhances  gas  exchange  (Zappa 
et  al.  2007),  and  generates  turbulence  that  mixes  the 
ocean  surface  layer  (Burchard  et  al.  2008;  Kukulka  and 
Brunner  2015).  Previous  observations  of  wave¬ 
breaking  turbulence  have  shown  strong  enhancement 
near  the  surface,  at  values  that  far  exceed  those  pre¬ 
dicted  by  simple  “law  of  the  wall”  boundary  layer 
scaling  (Agrawal  et  al.  1992;  Terray  et  al.  1996; 


Here,  the  wave  energy  loss  is  integrated  over  a  range  of 
frequencies  /,  and  the  vertical  reference  frame  for  depth 
z  can  either  be  wave  following,  noted  as  zw,  or  fixed, 
noted  as  z-  On  average,  the  wave-following  surface 
must,  by  definition,  become  the  mean  sea  level,  such  that 
(zw  =  0)  — ►  z  =  0.  However,  the  wave-resolved  s  may 
not  be  symmetric,  such  that  the  average  s  is  likely  dif¬ 
ferent  between  the  two  reference  frames.  Reconciling 
these  reference  frames  is  essential  to  comparing  the 


\Stodf~F~  je(z)dz 


(1) 


3  Denotes  Open  Access  content. 


Corresponding  author  address :  Jim  Thomson,  Applied  Physics 
Laboratory,  University  of  Washington,  1013  NE  40th  St.,  Seattle, 
WA  98105. 


E-mail:  jthomson@apl.washington.edu 


DOI:  10.1175/JPO-D-15-0130.1 


©  2016  American  Meteorological  Society 


1858 


JOURNAL  OF  PHYSICAL  OCEANOGRAPHY 


Volume  46 


different  observations  of  s  reported  in  the  literature,  and 
it  also  is  essential  to  obtaining  a  consistent  total  energy 
dissipation  rate  when  depth  integrating  \s(z)dz-  This 
simple  balance  neglects  transport  of  turbulence  and 
buoyancy  work  by  bubbles,  both  of  which  may  be  sig¬ 
nificant  in  strongly  forced  conditions. 

Formulations  for  the  total  energy  flux  F  typically  be¬ 
gin  with  the  momentum  flux  from  a  wind  stress  r  =  pau\a 
because  wave  breaking  in  deep  water  (i.e.,  whitecaps)  is 
dominated  by  the  short,  wind-driven  waves.  Here,  u*a  is 
the  air-side  friction  velocity,  which  we  convert  to  a 
water-side  friction  velocity  u*w  =  (pa/pw)ll2u*a,  using  the 
ratio  of  the  air  and  water  densities.  Dropping  the  sub¬ 
script  w ,  the  water-side  u *  is  used  throughout  the  for¬ 
mulations  that  follow.  The  wind  stress  acts  upon  the 
waves,  which  are  roughness  elements  at  the  free  surface, 
and  thus  some  formulations  for  total  energy  flux  include 
the  speeds  or  scales  of  the  waves  in  prescribing  the  work 
done  by  the  wind  stress.  The  subsequent  dissipation  of 
energy  happens  in  a  region  of  depth  z  where  waves  are  a 
dominant  process,  which  is  a  working  definition  of  the 
ocean  surface  layer  in  any  of  the  respective  vertical 
reference  frames. 

One  common  scaling  for  turbulence  in  the  ocean 
surface  layer  is  based  on  the  shear  production  of  an 
applied  stress  acting  on  a  surface  with  a  characteristic 
speed  (e.g.,  Craig  and  Banner  1994).  In  this  approach, 
the  input  of  turbulent  kinetic  energy  (TKE)  increases 
with  the  cube  of  the  friction  velocity 

F1  =  aul ,  (2) 

where  a  is  a  constant,  typically  O(102).  Although  the 
actual  values  exceed  the  log-layer  estimates,  the  physical 
scaling  for  the  TKE  input  is  a  generic  boundary  layer 
scaling  because  the  additional  velocity  scale  multiplied 
with  t  is  w*  (as  opposed  to  a  velocity  scale  associated 
with  the  waves).  Craig  and  Banner  (1994)  apply  this  u\ 
input  in  a  diffusive  model  for  the  vertical  structure  of  the 
turbulent  dissipation  rate  e(z),  which  is  now  in  regular 
usage  (e.g.,  Gerbi  et  al.  2013). 

Turbulence  generation  during  wave  breaking,  how¬ 
ever,  creates  a  boundary  layer  that  may  not  scale  as 
simply  as  the  cube  of  the  friction  velocity.  Wave 
breaking  is  an  intermittent  process  that  injects  turbu¬ 
lence  at  scales  related  to  the  waves  themselves,  rather 
than  the  wind  alone.  Recognizing  the  importance  of 
wave  scales,  Gemmrich  et  al.  (1994)  proposed  a  differ¬ 
ent  form  for  the  TKE  input,  which  uses  an  effective 
speed  ce  to  transfer  the  wind  stress  t  into  the  surface  as 
TKE,  such  that 


(3) 


Gemmrich  et  al.  (1994)  relate  ce  to  short  waves  that 
dominate  the  roughness  of  the  ocean  surface  and  find 
that  ce  ~  1ms-1.  This  is  a  bulk  value  that  is  meant  to 
represent  the  range  wave  speeds  that  carry  the  stress 
from  the  wind.  Although  the  correct  ce  is  still  an  open 
question,  the  important  aspect  of  this  formulation  is  that 
the  wind  stress  t  =  pau%  is  imparted  to  the  waves,  as 
represented  by  ce,  not  just  a  mean  sheared  current, 
which  would  be  represented  by  another  factor  of  u *  [as 
in  Eq.  (2)].  Conceptually,  ce  should  be  related  to  the 
speed  of  the  roughness  elements  (i.e.,  the  short  waves). 

Although  it  has  not  been  adopted  in  the  subsequent 
literature,  Phillips  (1985)  also  considered  an  adjust¬ 
ment  to  the  energy  fluxes  based  on  wave  scales.  Phillips 
(1985,  p.  522)  wrote,  “the  energy  source  for  near¬ 
surface  turbulence  represented  by  wave  breaking  is, 
in  mean,  distributed  over  a  wide  range  of  scales,  rather 
than  being  concentrated  at  the  energy-containing 
scales  as  it  is  in  shear-generated  turbulence.”  His  for¬ 
mulation  considers  the  available  flux,  defined  as  u%, 
that  can  be  carried  by  wavenumbers  bounding  the 
“equilibrium  range”  k0  and  k\.  The  result  is  a  formu¬ 
lation  for  a  balance  of  total  dissipation  and  total  TKE 
input  given  by 


F3  =  2yp3I(3p)ul  In 


(4) 


where  the  combined  constants  y(33I(3p)  are  approxi¬ 
mately  10-3.  The  wavenumber  limits  are  defined  by  a 
lower  bound  that  is  twice  the  peak  wavenumber  k0  =  2 kp 
and  an  upper  bound  that  is  caused  by  suppression  of 
waves  that  are  slower  than  the  drift  speed  of  the  surface, 
which  can  be  related  to  the  wind  friction  velocity  by 
k\  =  g/u%.  In  practice,  the  factor  ln^/Zco)  in  Eq.  (4) 
limits  the  increase  in  energy  flux  at  high  winds  by  nar¬ 
rowing  the  range  of  scales  on  which  the  wind  acts.  A 
wave  spectrum  with  a  smaller  equilibrium  range  will 
have  a  smaller  ratio  of  Aq/fco  and  thus  a  smaller  input  rate 
F3.  This  dependence  on  the  short  waves  is  loosely  similar 
to  using  an  effective  speed  that  remains  tied  to  short 
(and  relatively  slow)  waves  in  Eq.  (3). 

Phillips  (1985)  argued  that  these  wavenumbers  were 
the  scales  in  which  wave-breaking  dissipation  was  in 
balance  with  wind  input  and  nonlinear  transfers  within 
the  spectrum,  hence  the  term  equilibrium  range.  Equi¬ 
librium  also  was  proposed  by  Kitaigorodskii  (1983),  al¬ 
though  with  a  different  derivation.  The  equilibrium 
hypothesis  has  been  shown  to  be  consistent  with  wave¬ 
breaking  turbulence  measurements  under  moderate 
wind  speeds  (Thomson  et  al.  2013).  Assuming  equilib¬ 
rium,  and  using  a  wind  input  formulation  from  Plant 
(1982),  the  TKE  input  is 


F2  =  ceT  =  ceul. 
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^4  = 


=  J^cos(0)(27T f)E(f)  df,  (5) 


where  6  is  the  relative  difference  between  the  wind  and 
wave  direction  at  every  frequency  /  of  the  wave  energy 
spectrum  E ,  and  the  factor  of  0.04  in  Plant  (1982)  is 
omitted  because  the  adjustment  of  ( pJPw )1/2  has  already 
been  made  in  determining  the  water-side  u *.  In  deep 
water,  the  phase  speed  is  c  =  g/(27i/),  and  thus  the  net 
dependence  of  Swind  weights  the  high  frequencies  with 
/3.  This  is  similar  in  character  to  the  f4  weighting  of  the 
mean  square  slope  and  the  f  weighting  of  the  “wave 
saturation”  (Banner  et  al.  2002;  Hwang  et  al.  2012),  both 
of  which  have  been  correlated  with  wave-breaking  tur¬ 
bulence  (Gemmrich  2010). 

Finally,  dissipation  is  calculated  using  a  formulation  in 
the  WAVEWATCH  III  spectral  wave  model  (Tolman 
1991;  Tolman  and  Chalikov  1996;  Tolman  et  al.  2014). 
The  spectral  dissipation  from  the  wave  model,  in¬ 
tegrated  over  equilibrium  range  frequencies,  can  be 
used  as  yet  another  estimate  of  the  TKE  input: 


^5  = 


SbAdf, 


(6) 


which  is  applied  at  the  ocean  surface  wave  via  breaking. 

In  addition  to  the  total  energy  flux  and  equilibrium 
balance  of  breaking  waves,  the  depth  dependence  of  this 
balance  has  been  presented  in  several  previous  papers. 
For  example,  Terray  et  al.  (1996)  apply  the  model  for 
vertical  diffusion  of  the  TKE  input  from  Craig  and 
Banner  (1994)  and  assess  F1  and  F2  using  observations  of 
small  whitecaps  on  Lake  Ontario.  Terray  et  al.  (1996) 
evaluate  the  scaling 


(7) 


where  z  is  the  fixed  depth  beneath  the  mean  water  level,  A 
is  the  exponential  decay  rate,  and  Hs  is  the  significant  wave 
height.  Their  data  were  collected  in  conditions  of  7  to 
16  m  s_1  winds  and  0.2-  to  0.5-m  significant  wave  heights.  In 
terms  of  depth  dependence,  they  proposed  three  distinct 
vertical  regions  in  a  fixed  reference  frame:  1)  a  breaking 
zone,  0  <  z/Hs  <  0.6,  which  has  half  of  the  total  dissipation 
and  a  constant  profile;  2)  a  transition  layer,  which  follows 
(z/Hs)~2  and  depends  on  wave  age;  and  3)  a  deep  loga¬ 
rithmic  layer  ln(z).  They  did  not  make  measurements  in 
region  1  but  rather  inferred  this  using  the  dissipation 
necessary  to  meet  the  equilibrium  assumption. 

More  recently,  Gemmrich  (2010)  presented  mea¬ 
surements  in  a  wave-following  reference  frame  zw  and 
found  the  largest  dissipation  rates  to  occur  exclusively  in 


the  wave  crests,  above  the  still  water  level.  He  fit  the 
vertical  decay  of  the  turbulence  using  Eq.  (7)  and 
obtained  —1.5  <  A  <  —  1,  with  a  deep  log  layer  below 
one  Hs.  Similarly,  wave  tank  measurements  of  spilling 
waves  on  a  sloping  beach  suggest  that  80%  of  the  total 
dissipation  occurs  above  the  still  water  level  (Govender 
et  al.  2004).  Gemmrich  (2010)  converted  the  wave¬ 
following  measurements  to  a  fixed  vertical  coordinate, 
which  resulted  in  maximum  dissipation  at  z  ***  0.3 Hs. 
Since  the  adjusted  reference  frame  produces  an  appar¬ 
ent  decrease  in  the  near-surface  dissipation,  Gemmrich 
(2010)  noted  the  reference  frame  conversion  as  a  severe 
limitation  for  representing  wave-breaking  turbulence  in 
numerical  models  that  use  a  flat  surface.  The  Gemmrich 
(2010)  results  are  scaled  by  wind  stress  [Eq.  (2)]  in 
conditions  that  range  from  0  to  15  ms-1  winds,  with 
fetch-limited  waves  that  are  less  than  Hs  <  0.5  m.  Re¬ 
cently,  Sutherland  and  Melville  (2015)  used  similar 
methods  and  found  that  A  ~  1  near  the  surface  (within 
one  wave  height)  and  A  ~  2  below  that  [Eq.  (7)]. 

This  paper  uses  a  comprehensive  dataset  to  evaluate 
existing  formulations  for  wave-breaking  turbulence. 
Data  are  from  the  vicinity  of  Ocean  Weather  Station  P 
(OWS-P)  in  the  North  Pacific  (50°N,  145°W),  and  this 
work  extends  the  results  of  a  previous  study  at  that  site 
with  lower  winds  (Thomson  et  al.  2013).  The  previous 
study  analyzed  wave  equilibrium  but  did  not  scale  the 
surface  turbulence  nor  address  the  reference  frames. 
Here,  the  emphasis  is  scaling  the  surface  turbulence  in  a 
new  dataset  that  extends  to  20 ms-1  winds  and  recon¬ 
ciling  the  wave-following  reference  frame  of  the  obser¬ 
vations  with  the  fixed  reference  frame  of  the  models. 
Section  2  describes  the  data  collection  and  processing 
methods.  Section  3  presents  the  results.  Section  4  dis¬ 
cusses  the  implications.  Section  5  concludes. 


2.  Methods 

a.  Observations 

Wave,  wind,  and  surface  turbulence  data  were  col¬ 
lected  in  the  vicinity  of  OWS-P  at  50°N,  145°W.  The  data 
were  collected  during  a  mooring  turnaround  cruise  in 
January  2015  aboard  the  R/V  Thomas  G.  Thompson. 
Several  surface  wave  instrument  floats  with  tracking 
(SWIFTs;  see  Thomson  2012)  were  deployed  to  measure 
profiles  of  the  near-surface  turbulence  dissipation  rate. 
A  three-axis  sonic  anemometer  (RM  Young  model 
8100)  was  mounted  to  the  jackstaff  at  the  bow  of  the  ship 
to  measure  turbulent  winds  at  16-m  height  above  the 
surface.  Wave  spectral  measurements  were  collected 
from  the  SWIFTs  as  well  as  a  0.9-m  Datawell  directional 
waverider  (DWR  MKIII)  moored  at  OWS-P.  Video 
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Fig.  1.  Range  of  conditions  observed,  (a)  Wind  speeds,  (b)  significant  wave  heights,  and 
(c)  peak  wave  periods.  Black  points  are  from  SWIFTs,  green  crosses  are  from  the  shipboard 
sonic  anemometer,  and  blue  circles  are  from  the  waverider. 


measurements  of  wave  breaking  were  also  collected  and 
are  reported  in  a  separate  paper,  along  with  a  more 
detailed  analysis  of  the  wave  spectra  (Schwendeman  and 
Thomson  2015). 

Figure  1  shows  the  wind  and  wave  conditions  during 
the  cruise.  The  waves  were  predominantly  pure  seas, 
which  evolved  during  the  storms.  There  was  a  limited 
amount  of  swell  from  remote  sources.  Wave  conditions 
are  shown  from  both  the  SWIFTs  and  waverider;  how¬ 
ever,  the  waverider  values  are  only  shown  when  the 
SWIFTs  and  the  ship  were  within  30  nautical  miles  of  the 
waverider.  Using  a  10-min  burst  length  for  ensemble 
estimation,  the  full  dataset  is  composed  of  2522 
observations. 

b.  Wind  stress  estimates  u*a 

The  sonic  anemometer  data  from  the  jackstaff  of  the 
ship  were  processed  according  to  the  inertial  dissipation 
method  of  Yelland  et  al.  (1994),  in  which  an  air-side 
dissipation  rate  is  estimated  from  turbulence  spectra  and 
then  used  to  infer  the  wind  friction  velocity  u*a.  The 
sonic  anemometer  data  were  collected  at  10  Hz  and 
parsed  into  128-point  windows  that  were  despiked  fol¬ 
lowing  Goring  and  Nikora  (2002);  then  they  were  ta¬ 
pered  and  overlapped  50%  and  finally  fast  Fourier 
transformed.  Ensemble  spectra  were  made  at  10-min 
intervals  by  averaging  46  subwindows  to  obtain  final 


spectra  with  0.0391-Hz  frequency  resolution  and  92  de¬ 
grees  of  freedom. 

The  ensemble  spectra  were  fit  to  an  expected  fre¬ 
quency  dependence  of  /^5/3  in  the  inertial  frequency 
subrange  (1  <  /  <  4  Hz),  and  the  air-side  friction  velocity 
was  estimated  assuming  advection  of  a  frozen  field 
(Taylor’s  hypothesis)  at  a  speed  U,  such  that 


where  K  =  0.55  is  the  horizontal  Kolmogorov  constant, 
k  =  0.4  is  the  von  Karman  constant,  and  za  =  16  m  is  the 
measurement  height  above  the  still  water  level.  This 
assumes  neutral  stability,  which  is  justified  by  the  simi¬ 
larity  of  water  and  air  temperatures  during  the  obser¬ 
vations  (typically  within  2°C).  The  resulting  estimates  of 
wind  stress  r  =  pau%a  are  typically  within  10%  of  esti¬ 
mates  generated  using  standard  drag  coefficients  in 
t  =  paCD  U\ 0  (Smith  1988). 

The  wind  stress  estimates  from  the  ship  are  expanded 
using  the  wind  speed  measurements  from  the  SWIFTs 
and  standard  drag  coefficients,  which  are  applied  after 
adjusting  the  SWIFT  winds  from  za  =  lm  measurement 
height  to  the  =  10  m  standard.  The  air-side  friction 
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velocities  are  converted  to  water-side  friction  velocities 
w*  =  (pJPw)12u*cn  using  the  ratio  of  the  air  and  water 
densities. 

c.  WAVEWATCH  III  model  runs 

The  numerical  wave  model  used  here  is  WAVEWATCH 
III  (Tolman  1991;  Tolman  et  al.  2014),  with  the  physics 
package  of  Ardhuin  et  al.  (2010).  The  geographic  grid 
is  global,  at  0.5°  resolution,  and  the  spectral  grid  in¬ 
cludes  36  directional  bins  and  31  frequency  bins 
(0.0418  to  0.72  Hz,  logarithmically  spaced).  An  ob¬ 
struction  grid  is  used  to  represent  unresolved  islands 
(Tolman  et  al.  2014).  Winds  and  ice  concentrations  are 
taken  from  the  Navy  Global  Environmental  Model 
(NAVGEM;  Hogan  et  al.  2014),  with  the  former  input  at 
3-h  intervals  and  the  latter  at  12-h  intervals.  The  simula¬ 
tion  is  initialized  0000  UTC  1  December  2014  and  ends 
0000  UTC  1  February  2015.  The  required  spinup  time, 
during  which  predictions  of  remote  swells  are  invalid,  is 
estimated  to  end  0000  UTC  15  December  2014.  The 
physics  package  of  Ardhuin  et  al.  (2010)  requires  specifi¬ 
cation  of  a  parameter  /3max,  which  is  used  to  compensate 
for  the  mean  bias  of  the  input  wind  fields  or  lack  thereof; 
Pmax  =  1.2  is  used  for  this  hindcast.  The  model  version 
number  employed  is  development  version  5.08.  In  context 
of  the  present  application,  this  is  not  substantially  different 
from  public  release  version  4.18. 

d.  Wave  displacement  estimates 

The  SWIFTs  collected  GPS  velocities  of  4  Hz  on  a 
12-min  duty  cycle,  in  which  data  were  collected  for  a 
512-s  burst  and  then  processed  on  board  the  buoy  and 
transmitted  in  the  remaining  208  s  before  resuming  data 
collection.  Wave  energy  spectra  E(f)  were  estimated 
following  the  methods  of  Herbers  et  al.  (2012).  The 
waverider  collected  buoy  pitch,  roll,  and  heave  dis¬ 
placements  at  1.28  Hz  on  a  30-min  duty  cycle,  following 
the  Datawell  standards.  Both  types  of  buoys  processed 
data  on  board  at  the  end  of  each  duty  cycle  and  trans¬ 
mitted  the  spectral  moments  to  servers  on  shore  via 
Iridium  satellites. 

In  addition  to  the  GPS  receivers,  the  SWIFT  buoys 
were  equipped  with  MicroStrain  3DM-GX3  -35  sensors, 
which  measure  three  axes  of  acceleration,  rotation  rate, 
and  orientation  at  25  Hz.  These  measurements  were 
used  to  define  the  quaternion  orientation  matrix  at 
25  Hz,  from  which  true  vertical  heave  accelerations  (in 
the  earth  reference  frame)  were  calculated.  These  ac¬ 
celerations  were  integrated  twice  in  the  time  domain, 
with  a  high-pass  filter  at  each  integration,  to  reconstruct 
a  wave-resolved  time  series  of  sea  surface  elevations  77, 
relative  to  mean  sea  level.  The  elevations  were  used  to 
calculate  wave  energy  spectra,  and  these  are  compared 


with  the  energy  spectra  from  the  GPS  velocity  method 
of  Herbers  et  al.  (2012)  and  with  the  output  of  the 
Datawell  waverider.  Agreement  is  within  5%  for  most 
cases.  These  elevations  were  down  sampled  to  0.5  Hz 
and  used  to  map  turbulence  measurements,  collected  in 
the  wave-following  reference  frame  of  the  buoy  zw,  to  a 
fixed  reference  frame  z- 


e.  Turbulent  dissipation  rate  estimates  s 

The  SWIFTs  collected  pulse-coherent  Doppler  sonar 
(Nortek  Aquadopp  HR)  profiles  of  turbulent  velocities 
beneath  the  wave-following  surface  denoted  as  =  0. 
The  details  of  data  collection,  quality  control,  and  pro¬ 
cessing  are  described  in  Thomson  (2012)  and  will  only 
be  reviewed  here.  The  primary  differences  here  are  to 
use  the  double-sided  structure  function  (i.e.,  velocity 
variations  above  and  below  each  bin)  and  to  calculate 
wave-resolved  estimates  before  creating  burst-averaged 
ensemble  values.  The  wave-resolved  estimates  are  nec¬ 
essary  to  directly  map  the  measurements  from  a  wave¬ 
following  reference  frame  to  a  fixed  reference  frame. 

The  turbulent  velocities  were  collected  at  4  Hz  and 
were  processed  to  estimate  the  second-order  structure 
function  of  the  turbulent  velocity  fluctuations  beneath 
the  wave-following  surface  at  2-s  intervals  (=0.5  Hz) 
after  application  of  pulse  correlation  quality  control 
(minimum  correlation  of  30).  The  structure  function  is  a 
direct  spatial  realization  of  the  theoretical  Kolmogorov 
(1941)  energy  cascade  from  large  to  small  scales.  Fitting 
the  observed  structure  function  D(zw,r)  ~  r2/3,  where  r  is 
the  spatial  separation  of  velocity  measurements  along  a 
profile,  is  equivalent  to  fitting  a  k~5/ 3  wavenumber 
spectrum,  and  thus  the  turbulent  dissipation  rates  were 
estimated  according  to  Wiles  et  al.  (2006): 


®(0 = 


AVJ 


(9) 


where  Cv  =  1.45  is  a  constant,  and  A(zw)  is  the  ampli¬ 
tude  of  the  structure  function  determined  for  each 
from  fitting  D(zw,r)  =  Ar2/ 3  -  N.  Here,  N  is  the  con¬ 
tribution  to  variance  from  noise  in  the  Doppler  mea¬ 
surements;  it  is  a  free  parameter  in  the  fit  to  r2/3  and  is 
checked  against  the  expected  noise  of  the  profiler  (see 
Thomson  2012). 

f  Vertical  reference  frames 

The  wave-resolved  values  (2  s)  are  mapped  to  the 
fixed  reference  frame  and  used  to  calculate  burst- 
averaged  values  (512  s).  In  the  wave-following  frame, 
=  0  is  the  instantaneous  surface,  and  results  go  down 
to  zw  =  —0.5  m.  In  the  fixed  reference  frame,  z  =  0  is  the 
mean  sea  level,  and  77  is  the  instantaneous  wave  surface 
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Fig.  2.  (a)-(d)  An  example  30-s  time  series  of  a  breaking  wave,  which  is  confirmed  using 
images  collected  on  board  the  SWIFT.  The  wave-resolved  turbulent  dissipation  rate  profiles 
are  measured  in  a  (e)  wave-following  reference  frame  that  is  directly  mapped  to  a  (g)  fixed 
reference  frame  using  the  surface  elevation  time  series,  burst-averaged  profiles,  using  512  s  of 
data,  and  (f),(h)  shown  for  each  reference  frame.  The  fixed  reference  frame  results  include  data 
extrapolation  from  the  wave-following  reference  frame,  such  that  all  depths  are  averaged 
uniformly.  The  dashed  line  in  (f)  shows  a  burst-averaged  dissipation  rate  profile  calculated 
directly  from  all  the  velocities  measured  in  the  burst  rather  than  taking  the  average  of  the  wave- 
resolved  profiles. 


relative  to  mean  sea  level.  The  wave-resolved  e(zw) 
values  are  mapped  to  s(z)  using  z  =  v  -  zw,  where  z  is 
gridded  at  0.1 -m  resolution.  Both  e(zw)  and  s(z)  are 
burst  averaged  over  512  s  at  each  value  of  zw  or  z  to 
create  (e(zw))  and  (s(z)). 

In  mapping  to  the  gridded  z  reference  frame,  only  a 
subset  of  the  grid  cells  is  populated.  To  avoid  biases 
caused  by  data  sparsity  when  subsequently  calculating 
burst-averaged  profiles  of  dissipation  in  the  fixed  refer¬ 
ence  frame,  the  missing  values  are  filled  via  extrapola¬ 
tion.  The  extrapolation  is  done  in  the  wave-following 
reference  frame  before  mapping  to  the  fixed  reference 
frame.  The  wave-resolved  dissipation  values  above  sea 
surface  are  assigned  zero.  The  dissipation  values  below 
the  deepest  observation  are  extrapolated  by  fitting  the 
power-law  decay  of  Eq.  (7)  using  a  constant  A  =  1.4, 
which  is  determined  from  a  logarithmic  fit  to  the  entire 
dataset. 

Figure  2  shows  an  example  of  the  wave-resolved  dis¬ 
sipation  rates  and  the  corresponding  burst-averaged 
results  in  each  reference  frame.  This  particular  exam¬ 
ple  has  a  breaking  wave  at  t  =  113  s,  as  indicated  by  the 
surface  images  collected  on  board  the  buoy  every  4  s 


(Figs.  2a-d).  Note  that  only  30  s  of  data  are  shown  in  the 
wave-resolved  panels,  such  that  the  individual  breaker 
can  be  seen,  but  the  burst-averaged  panels  show  results 
calculated  using  the  whole  512-s  data  burst.  Other  bursts 
are  similar,  including  the  shift  of  the  very  shallow  s(zw) 
profiles  to  much  greater  depths  in  s(z).  The  shift  occurs 
because  the  changing  surface  elevation  77  carries  the  zw 
coordinates  through  several  meters  of  z.  Another  simi¬ 
larity  across  the  various  bursts  of  data  is  the  long  time 
scale  of  e,  which  remains  elevated  after  a  breaking  wave 
for  more  than  one  wave  period.  This  means  that  the 
strong  turbulence  generated  during  breaking  has  time  to 
be  carried  through  a  wide  range  of  fixed  z  values, 
nominally  —HJ2  <  z  <  HJ2 ,  before  returning  to  a 
background  level.  The  effect  is  to  spread  the  high 
e  values  through  the  surface  layer  rather  than  concen¬ 
trate  them  exclusively  above  mean  sea  level. 

The  maximum  dissipation  in  the  fixed  reference  frame 
occurs  near  the  mean  sea  level  (z  =  0)  and  decreases 
both  above  and  below  this  level  (Fig.  2h).  Other  bursts 
are  similar,  and  the  general  result  is  that  the  average 
turbulence  accumulates  at  the  average  surface.  Again, 
this  is  because  the  turbulence  persists  longer  than  wave 
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Fig.  3.  Time  series  of  processed  results  (a)  wind  friction  velocity,  (b)  burst-averaged  tur¬ 
bulent  dissipation  rate  in  the  wave-following  reference  frame,  and  (c)  burst-averaged  turbulent 
dissipation  rate  in  the  fixed  reference  frame. 


phase  changes.  Above  mean  sea  level,  the  decrease  in 
the  average  is  caused  by  assigning  zero  dissipation  when 
the  grid  does  not  have  water.  Below  mean  sea  level,  the 
decrease  in  the  average  is  caused  by  the  observed  and 
extrapolated  decay  of  turbulence  [Eq.  (7)]. 

3.  Results 

The  processed  wind  and  turbulence  results  are  shown 
in  Fig.  3  for  a  total  of  2522  bursts  of  data,  each  a  burst- 
averaged  ensemble  lasting  512  s.  Many  of  these  bursts 
overlap  in  time  because  multiple  SWIFTs  were  de¬ 
ployed  and  collected  data  simultaneously.  The  wind 
friction  velocity  estimates  from  the  sonic  anemometer 
are  consistent  with  application  of  a  drag  coefficient  to 
the  SWIFT  wind  speed  measurements;  both  increase 
monotonically  to  approximately  u*a  «  1  ms-1  for  the 
max  observed  wind  speed  of  £/i0  =  20  m  s-1.  The  air-side 
value  u*a  presented  in  this  figure  is  converted  to  a  water¬ 
side  value  u *  =  ( pjpw)ll2u*a  for  the  remaining  pre¬ 
sentation  of  results  (i.e.,  this  figure  is  the  only  place 
where  u*a  appears). 

In  Figs.  3b  and  3c,  the  burst-averaged  turbulent  dis¬ 
sipation  rates  in  the  wave -following  and  fixed  reference 
frames  track  the  time  series  of  wind  forcing.  The  values 
are  generally  in  the  range  of  10-4  <  s  <  l(T2ni2s“3, 


although  some  values  are  as  low  as  10-5m2s-3  during 
the  very  calm  period  on  10  January  2015.  The  range  is 
consistent  with  previous  observations  of  the  turbulent 
dissipation  rate  in  the  surface  layer  (e.g.,  Gemmrich 
2010;  Schwendeman  et  al.  2014;  Sutherland  and  Melville 
2015).  In  the  wave-following  reference  frame  of  Fig.  3b, 
the  dissipation  rates  are  maximum  at  the  surface  {zw  ~  0) 
and  decrease  with  depth.  In  the  fixed  reference  frame  of 
Fig.  3c,  the  dissipation  rates  are  maximum  at  the  mean 
sea  level  (z  =  0)  and  decrease  over  several  meters  in 
depth.  This  is  because  the  time  scales  for  turbulent 
dissipation  are  longer  than  the  time  scales  of  individual 
waves.  As  shown  in  the  example  of  Fig.  2,  the  turbu¬ 
lence  associated  with  an  individual  breaking  wave  per¬ 
sists  for  more  than  one  wave  period.  This  means  the 
turbulence  can  be  carried  from  the  crest  level  down  to 
the  trough  level  in  z,  even  though  it  may  be  isolated  very 
near  the  surface  in  zw-  Thus,  the  time  series  of  e(z)  in 
Fig.  3c  populate  a  larger  range  of  the  z  grid  when  the 
waves  are  larger. 

Figure  3  has  scatter  in  the  dissipation  rates.  The  values 
at  the  highest  wind  forcing  are  mixed  with  some  of  the 
moderate  conditions.  Some  of  this  can  be  attributed  to 
the  episodic  nature  of  breaking  and  probability  of  the 
number  and  size  of  individual  breaking  waves  within 
each  ensemble  (a  burst-averaged  result  from  512  s)  used 
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Dissipation  rate,  f  (e(zw))dzw  —  f  e^dz 


Fig.  4.  Integrated  rate  of  energy  input  from  the  wind  using  five  different  formulations  vs 
measured  total  turbulence  dissipation  rate  in  the  upper-ocean  surface  layer.  Symbols  show  the 
bin  averages  and  vertical  lines  show  the  standard  errors. 


for  SWIFT  calculations.  Of  course,  wave  dynamics,  in 
addition  to  wind  forcing,  control  wave  breaking,  and  the 
combined  effects,  are  evaluated  using  the  various  for¬ 
mulations  for  TKE  flux  into  the  ocean  surface. 

a.  Integrated  (total)  energy  fluxes 

The  assumption  in  most  treatments  of  wave-breaking 
turbulence  is  that  the  input  rate  of  TKE  at  the  ocean 
surface  matches  the  dissipation  rate  below  the  surface 
[Eq.  (1)].  Figure  4  compares  the  five  different  estimates 
for  the  total  TKE  input  rate  F  with  the  observed  tur¬ 
bulence  dissipation  rates  [depth-integrated  total  in  the 
wave-following  reference  frame  without  extrapolation; 
J (e(zw))  dzw]  using  binned  averages.  The  comparisons 
are  summarized  in  Table  1,  which  shows  the  raw  corre¬ 
lation  coefficients  from  all  2522  ensembles  (rather  than 
the  bin  averages).  Here,  we  use  a  constant  a  =  75  in  F1 
[Eq.  (2)]  and  a  constant  transfer  velocity  of  ce  =  2 ms-1 
in  F2  [Eq.  (3)].  These  were  selected  as  the  values  within 
the  range  from  the  literature  for  young  seas  that  give 
the  best  comparison  with  the  observations  of  dissi¬ 
pation.  The  WAVEWATCH  III  dissipation  values 
were  obtained  by  integrating  §Sbrkdf  over  frequencies 
greater  than  twice  the  peak  frequency.  This  frequency 


integration  limit  is  chosen  because  swell  frequencies 
have  a  different  dissipation  mechanism  (Ardhuin  et  al. 
2010). 

For  the  comparison,  a  background  J  (e0)  dzw  =  0.7  X 
10-3m3s-3  is  removed  from  each  ensemble.  This  is 
consistent  with  previous  studies  (e.g.,  Thomson  et  al. 
2009;  Gemmrich  2010;  Schwendeman  et  al.  2014)  where 
the  intent  is  to  isolate  the  wave-breaking  dissipation 
from  dissipation  due  to  shear-driven  turbulence 
(McWilliams  et  al.  2012),  microbreaking  (Sutherland 
and  Melville  2015),  and  swell  dissipation  (Ardhuin  et  al. 
2010).  Here,  we  have  used  the  wave-resolved  estimates 


Table  1.  TKE  flux  estimates  F  and  correlation  to  observed  dissipation 

J'e(-^vv)  dzw  • 


Flux  estimate 

Correlation  coefficient  R2 

Fi  =  OLU% 

0.48 

F2  =  ceu*2 

0.53 

F3  =  2y(l3I(3p)ui  \n(ki/ko) 

0.48 

r  u2 

F4  =  g^cos(0)(,27Tf)E(f)df 

0.49 

F5  =  J  Sbrk  df 

0.52 
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to  dissipation  rate  to  further  constrain  the  determination 
of  the  background  value.  We  define  background  as  the 
median  of  the  lowest  half  of  the  values  in  a  given  burst 
(512  s).  Note  that  the  dissipation  values  are  highly  non- 
Gaussian  because  very  large  values  occur  episodically 
during  breaking;  thus,  this  definition  is  not  equivalent  to 
the  median  of  the  values  less  than  the  mean.  Using  this 
definition,  the  background  dissipation  has  almost  no 
variation  as  a  function  of  wind  speed  above  5ms”1.  The 
background  values  are  in  the  range  of  0.5  X  1(T3  to  0.9  X 
10-3m3s-3,  with  minimal  dependence  on  wave  age  or 
wind  speed.  Thus,  we  use  a  constant  0.7  X  10-3  m3  s-3 
throughout  this  analysis.  Although  the  definition  using 
the  median  lowest  half  of  the  values  may  seem  an  arbi¬ 
trary  choice,  this  background  level  is  confirmed  as  the 
threshold  when  visible  breaking  becomes  the  dominant 
process,  following  Gemmrich  (2010). 

All  of  the  total  flux  estimates  show  general  agreement 
with  the  bin-averaged  measurements  of  turbulent  dis¬ 
sipation  rate.  The  formulation  using  a  transfer  velocity 
[F2  from  Eq.  (3)]  has  the  best  correlation,  although  the 
differences  among  the  correlation  coefficients  are  only 
marginally  significant  (80%  level).  The  general  balance 
of  input  and  dissipation  is  consistent  with  the  equilib¬ 
rium  concept  of  Phillips  (1985)  and  the  observations  of 
Thomson  et  al.  (2013),  although  that  study  was  at  lower 
wind  speeds.  Similarly,  the  persistence  of  a  /^4  slope  in 
the  wave  spectra  shown  in  Schwendeman  and  Thomson 
(2015)  also  suggests  equilibrium,  following  Babanin  and 
Soloviev  (1998). 

Relative  to  the  in  situ  dissipation,  most  of  the  esti¬ 
mates  of  the  TKE  input  in  Fig.  4  are  low  during 
moderate  conditions  (<0.7  X  10-3m3s-3)  and  high 
during  rough  conditions  (>0.7  X  10-3m3s-3).  The 
WAVEWATCH  III  estimates  of  §Sbrkdf  are  particu¬ 
larly  low  during  the  moderate  conditions.  This  might  be 
related  to  wave  growth  during  those  conditions,  although 
the  bulk  growth  rate  is  typically  less  than  10-4m2s-3. 
The  transfer  velocity  formulation  F2  [Eq.  (3)]  is  the  best 
for  matching  the  in  situ  dissipation  during  moderate 
conditions,  but  it  too  exceeds  the  in  situ  dissipation 
values  during  rough  conditions.  That  all  estimates  exceed 
the  measurements  during  rough  conditions  suggests  a 
process  is  missing  from  the  data  analysis  or  a  deficiency  in 
the  observational  data  itself.  It  may  be  that  these  data 
simply  do  not  extend  deep  enough  to  capture  all  of  the 
dissipation.  Alternatively,  bubbles  may  be  important, 
both  to  the  process  of  energy  dissipation  and  to  the 
quality  of  acoustic  Doppler  data,  especially  during  rough 
conditions. 

Previous  results  have  shown  that  the  work  done 
against  the  buoyancy  of  bubbles  in  breaking  waves  is 
additional  sink  of  energy  (e.g.,  Lamarre  and  Melville 


1991).  Furthermore,  in  a  competing  mechanism,  the 
presence  of  large  bubble  clouds  near  the  surface  may 
actually  limit  surface  turbulence  in  high  winds.  Deane 
et  al.  (2016)  suggest  that  TKE  dissipation  rates  are 
limited  to  0(10°)  m2s-3  by  the  buoyancy  stabilization 
within  bubbles  clouds  and  that  increasingly  energetic 
wave  breaking  simply  creates  larger  bubble  clouds  that 
achieve  this  limit.  The  highest  dissipation  rates  calcu¬ 
lated  here  do  not  approach  the  Dean  et  al.  values  but 
that  again  may  be  a  limitation  of  the  Doppler  sonars. 

The  relationship  of  bubbles  and  turbulence  is  also 
highlighted  by  the  recent  work  of  Lim  et  al.  (2015),  who 
show  that  accounting  for  void  fraction  is  crucial  to  in¬ 
terpreting  the  rapid  dissipation  of  wave  energy 
beneath  a  breaking  wave.  For  the  roughest  conditions  in 
this  dataset  (20 ms-1  winds),  bubbles  and  “spindrift” 
(spraying  foam)  may  become  important  to  the  total 
energy  flux  budget.  For  spindrift  in  particular,  the 
transfer  of  kinetic  energy  effectively  bypasses  the  waves 
because  energy  does  not  get  input  to  the  wave  spectrum 
(where  it  would  be  available  to  dissipate  in  the  breaking 
process).  Restated,  the  “clipping”  off  the  top  of  wave 
crests  is  a  direct  transfer  of  energy  on  a  wave-by-wave 
basis,  as  opposed  to  the  indirect  transfer  described  by 
the  equilibrium  of  wind  input  and  breaking  dissipation 
achieved  via  the  exchanges  of  many  waves  in  a  spectrum. 

Bubble  dynamics  are  beyond  the  scope  of  this  paper, 
and  spray  is  only  beginning  to  occur  at  the  upper  end 
of  this  dataset  (U10  =  20 ms-1).  The  pulse-coherent 
acoustic  Doppler  methods  used  on  board  the  SWIFTs 
are  not  capable  of  measuring  turbulent  velocities  inside 
bubble  clouds.  The  bubble  clouds  are  screened  during 
initial  quality  control  of  the  raw  data,  using  the  pulse 
correlations  and  acoustic  backscatter,  which  are  reduced 
and  elevated,  respectively,  by  bubbles.  Thus,  the  results 
herein  are  confined  to  the  average  TKE  dissipation  rates 
outside  of  the  bubble  clouds.  This  might  significantly 
change  interpretation  of  the  total  energy  flux  balance  at 
the  surface.  However,  the  active  bubble  clouds  cover 
less  than  1%  of  the  sea  surface  at  any  given  time,  even  at 
winds  of  >15  m  s-1  (Schwendeman  and  Thomson  2015), 
and  thus  our  estimates  of  the  burst-averaged  TKE  dis¬ 
sipation  rate  (e)  can  be  interpreted  as  values  that  are 
characteristic  for  most  of  the  ocean  surface  layer.  Even 
though  the  total  dissipation  obtained  by  these  methods 
may  be  biased  low,  these  characteristic  values  may  be 
more  appropriate  when  using  dissipation  as  a  proxy  for 
gas  transfer  (Zappa  et  al.  2007). 

b.  Depth  scaling 

In  the  wave-following  reference  frame,  the  strongest 
turbulent  dissipation  rates  are  isolated  to  a  region  very 
close  to  the  surface  (zw  >  —0.2  m),  even  when  the  waves 
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are  several  meters  in  height.  This  is  consistent  with 
several  published  observations  (Gemmrich  2010; 
Thomson  2012;  Sutherland  and  Melville  2015).  How¬ 
ever,  other  published  observations,  using  a  fixed  refer¬ 
ence  frame,  have  concluded  that  strong  dissipation 
occurs  much  farther  below  the  surface  (Terray  et  al. 
1996;  Drennan  et  al.  1996;  Feddersen  2012).  The  ap¬ 
parent  contradiction  can  be  reconciled  by  considering 
the  effect  of  wave  orbital  motions  in  converting  the 
wave-following  reference  frame  to  a  fixed  reference 
frame  using  wave-resolved  observations.  The  orbital 
motions  move  the  surface  (z  =  rj,  zw  =  0)  from  crest  to 
trough  and  back  every  wave  period.  While  zw  =  0  re¬ 
mains  the  same,  the  fixed  reference  frame  z  ranges  from 
7]  ~  —HSI2  to  7]  ***  HJ2.  If  the  strong  turbulence  of  a 
breaking  crest  persists  more  than  half  a  wave  period,  it 
will  be  carried  down  to  the  trough  level  z  =  —HSI2  while 
still  appearing  to  be  isolated  within  a  shallow  region 
near  the  surface,  -0.4  <  zw  <  0  m,  as  shown  in  the  ex¬ 
ample  of  Fig.  2. 

Here,  the  wave-resolved  (at  2-s  time  steps)  estimates 
of  dissipation  rate  are  converted  using  the  direct  trans¬ 
form  z  =  7]  -  zw  and  then  burst  averaged  (at  512-s  time 
steps).  The  burst-averaged  values  are  denoted  with 
brackets.  The  burst-averaged  profiles  in  both  the  zw  and 
z  reference  frames  are  then  used  to  evaluate  the  depth 
scaling  of  Eq.  (7),  which  uses  significant  wave  height  to 
normalize  a  power  law  decreasing  in  depth.  A  single 
formulation  of  TKE  flux  F2  is  used,  which  is  selected 
based  on  the  quality  of  correlation  in  the  total  dissipa¬ 
tion  comparison  (Table  1).  Using  other  flux  estimates 
has  a  negligible  effect  on  the  resulting  depth  de¬ 
pendence  because  the  flux  estimates  are  similar  over 
most  of  the  range  of  the  observations. 

The  depth-scaled  results  are  presented  in  Fig.  5  and 
Table  2.  In  keeping  with  the  original  formulation,  the 
scaled  results  are  only  presented  for  depths  below  mean 
sea  level,  and  it  is  important  to  note  that  this  misses 
approximately  half  of  the  total  dissipation.  The  wave¬ 
following  zw  results  in  decay  below  the  surface  and  fit  a 
power  law  with  an  exponent  of  approximately  A  ^  1.4. 
This  is  consistent  with  other  wave-following  observa¬ 
tions  (Gemmrich  2010;  Zippel  and  Thomson  2015; 
Thomson  et  al.  2014)  and  region  2  of  the  Terray  et  al. 
(1996)  observations  in  a  fixed  reference  frame.  The  fixed 
z  results  have  a  similar  depth  dependence,  though  shif¬ 
ted  deeper,  and  the  exponent  is  A  «  1.6.  Both  reference 
frames  suggest  a  region  of  quasi-constant  dissipation 
near  the  surface,  which  was  postulated,  but  not  mea¬ 
sured,  in  the  fixed  frame  (Terray  et  al.  1996,  their  region 
1).  In  the  wave -following  frame,  the  constant  region  is 
only  1%  of  the  normalized  depth  (zJHs  >  -1CT2).  In 
the  fixed  reference  frame,  the  constant  region  is  at  least 


Fig.  5.  Scaled  vertical  profiles  of  turbulent  dissipation  rates  in 
(a)  the  wave-following  reference  frame  and  (b)  the  fixed  reference 
frame.  Diamonds  show  bin  averages,  and  horizontal  lines  show  the 
standard  percent  errors  (which  display  symmetrically  in  logarithm 
space)  at  each  bin. 
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Table  2.  Timescales,  reference  frames,  and  scaled  results  for  the  measured  profiles  of  TKE  dissipation  rate  s.  Confidence  intervals  on  the 

scaling  fits  are  at  95%  significance  level. 


Dissipation 

estimate 

Time  scale 

Reference  frame 

Scaling 

b(zw) 

s(z) 


Phase  resolved  (2  s) 

Wave 

Phase  resolved  (2  s) 

Fixed 

Burst  averaged  (512  s) 

Wave 

following  (zw  =  0  at  sea  surface) 
(z  =  0  at  mean  sea  level) 

following  (zw  =  0  at  sea  surface) 


N/A 

N/A 


— 1.41±0.02 


<10~2 


(s(z)> 


Burst  averaged  (512  s)  (z  =  0  at  mean  sea  level) 


— 1.56±0.03 

for  -4-<10‘1 

Ms 


10%  of  the  normalized  depth  (z/Hs  >  -1CT1).  The 
measurements  of  the  present  study  do  not  extend  deep 
enough  to  observe  the  third,  and  deepest,  layer  dis¬ 
cussed  in  Terray  et  al.  (1996),  which  is  a  log  layer  be¬ 
neath  the  wave-affected  layer. 

Although  these  results  reconcile  much  of  the  differ¬ 
ences  in  the  literature  between  fixed  reference  frames 
(e.g.,  Terray  et  al.  1996)  and  wave-following  reference 
frames  (e.g.,  Gemmrich  2010),  important  distinctions 
remain.  Terray  et  al.  (1996)  extrapolated  dissipations 
above  the  trough  level  [zJHs  >  O(-10°)]  such  that  the 
total  depth-integrated  dissipation  matched  the  TKE 
input  rate  F  and  the  dissipation  above  the  mean  sea  level 
was  zero.  As  there  are  two  parameters  (e  and  z)  and  only 
one  hard  constraint  (matching  F),  there  are  many  other 
extrapolations  that  would  match  F  and  yet  allow  non¬ 
zero  dissipation  above  the  mean  sea  level.  Revisiting 
Fig.  7  of  Terray  et  al.  (1996),  it  is  clear  that  the  extrap¬ 
olation  is  not  well  constrained  because  there  are  only 
four  observations  above  [zw/Hs  >  0(  — 10°)].  Below 
their  extrapolation,  the  observations  of  Terray  et  al. 
(1996)  are  roughly  consistent  with  the  observations  of 
the  present  study. 

4.  Discussion 

In  summary,  the  enhanced  turbulence  dissipation 
rates  associated  with  wave  breaking  are  isolated  in  an 
extremely  shallow  surface  layer  of  h  <  1  m,  yet  they  also 
are  distributed  throughout  a  thicker,  wave-affected 
layer  of  h  <  Hs.  The  answer  depends  on  the  reference 
frame.  For  the  purpose  of  comparing  with  numerical 
circulation  models,  which  do  not  typically  include  wave- 
resolved  surface  waves,  the  fixed  reference  frame  is 
commonly  used  (e.g.,  Gerbi  et  al.  2013)  and  the  present 
results  indicate  that  turbulent  dissipation  should  be  in¬ 
cluded  for  depths  h  ~  Hs.  The  lingering  problem  is  then 
what  to  do  with  the  portion  of  total  dissipation  that  oc¬ 
curs  above  mean  sea  level  (z  >  0).  For  this  dataset,  once 
converted  to  a  fixed  reference  frame,  that  portion  is 
approximately  half  of  the  total. 


The  observed  quasi-constant  layer  of  dissipation 
within  the  upper  10%  of  the  nondimensional  fixed 
reference  frame  ( ~z/Hs  >  10_1)  is  qualitatively  con¬ 
sistent  with  the  recent  large-eddy  simulation  (LES) 
results  of  Kukulka  and  Brunner  (2015),  and  the  near- 
complete  dissipation  of  most  wave-breaking  energy 
within  the  upper  few  meters  is  consistent  with  Noh 
et  al.  (2004).  Results  are  also  consistent  with  a  diffusive 
layer  below  the  constant  layer  (e.g.,  Craig  and  Banner 
1994;  Terray  et  al.  1996).  Results  do  not  extend  deep 
enough  to  evaluate  the  expected  log-layer  below 
~z/Hs  <  0(10°).  The  time  scales  of  the  observed 
dissipation  are  also  qualitatively  consistent  with 
numerical  studies.  Specifically,  Sullivan  et  al.  (2004) 
used  direct  numerical  simulation  (DNS)  to  show  that  a 
single  breaking  wave  can  energize  the  surface  layer  for 
more  than  50  wave  periods,  and  this  long  time  scale  is 
central  to  the  vertical  advection  of  turbulence  by  orbital 
motions  and  phase-resolving  change  of  reference  frame 
presented  here. 

As  to  which  input  TKE  flux  to  use  in  numerical 
models,  the  Gemmrich  et  al.  (1994)  effective  transfer 
speed  formulation  [F2\  Eq.  (3)]  has  the  best  agreement 
with  these  observations,  although  only  marginally  so. 
All  formulations  exceed  the  observed  dissipation  at 
high  winds,  and  this  highlights  the  cubic  u%  de¬ 
pendence  of  F1  [Eq.  (2)].  The  Phillips  (1985)  formu¬ 
lation  [F3;  Eq.  (4)]  may  be  an  improvement  over  the 
simpler  u\,  but  it  is  important  to  note  that  F3  here  is 
really  just  u\  with  a  different  coefficient,  since  the  ratio 
h^/q/ko)  is  fixed  at  ln(2.5)  with  our  definition  of 
the  equilibrium  range  as  k0  <  km  <  k\.  If,  instead,  the 
k\  =  gu*2  prescription  from  Phillips  (1985)  is  used,  the 
total  TKE  input  rate  does  indeed  have  a  dependence 
less  steep  than  u%.  This  alternate  form  has  a  negligible 
effect  in  comparing  to  this  dataset;  however,  the  sub¬ 
tleties  of  the  Phillips  (1985)  formulation  may  be  useful 
in  future  studies. 

Assuming  a  total  energy  balance  [Eq.  (1)]  for  every 
burst-averaged  ensemble  (2522  total),  effective  transfer 
speeds  can  be  estimated  rather  than  prescribed  as  the 
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Fig.  6.  Empirical  transfer  velocities  ce  as  (a),(b)  histograms  and  (c),(d)  correlated  with  wave 
age.  Values  are  (left)  dimensional  and  (right)  scaled  by  peak  wave  phase  speed.  Young  waves 
(wave  age  <  2)  are  shown  in  blue.  Mature  seas  (wave  age  >  2)  are  shown  in  green.  The  constant 
ce  =  2  m  s_1  used  in  the  earlier  sections  is  shown  as  a  red  line  in  (a). 


constant  ce  =  2 ms-1  used  in  the  preceding  sections. 
Figure  6  shows  these  empirical  transfer  speeds  in  his¬ 
tograms,  and  versus  wave  age,  as  dimensional  and  non- 
dimensional  (normalized  by  peak  wave  phase  speed) 
values.  The  dimensional  histograms  are  consistent  with 
the  2 ms-1  value  taken  as  a  constant  in  the  preceding 
analysis.  This  also  is  consistent  with  the  concept  of  short 
waves,  which  travel  slowly,  providing  the  roughness  el¬ 
ements  upon  which  the  wind  stress  can  do  work  (i.e., 
transfer  energy). 

The  values  in  Fig.  6  are  separated  into  young  seas 
(blue  points)  and  mature  seas  (green  points)  using  the 
wave  age  (ratio  of  the  peak  wave  phase  speed  cp  to  the 
wind  speed  U10).  For  young  seas,  the  equilibrium  input- 
dissipation  balance  works  well,  and  the  histogram  is 
centered  on  ce  =  2 ms-1  as  expected.  When  normalized 
by  the  peak  phase  speed,  the  histogram  is  centered  on 
celcp  =  1/5.  For  mature  seas,  the  equilibrium  balance  is 
not  a  valid  assumption  because  dissipation  is  likely  to 
exceed  input.  Figure  6  shows  the  spuriously  high  ce 
values  that  result  in  order  to  match  large  dissipations 
with  smaller  inputs  if  the  equilibrium  calculation  is  ap¬ 
plied  to  mature  seas.  The  effect  is  most  severe  at  large 
wave  ages.  The  clear  recommendation  is  to  restrict 


application  of  the  effective  transfer  velocity  formulation 
to  young  seas. 

5.  Conclusions 

In  situ  observations  of  turbulent  dissipation  rates 
immediately  below  the  wave-following  surface  of  the 
ocean  are  consistent  with  several  estimates  for  the  input 
of  TKE  by  a  wind  stress  performing  work  on  a  wave 
surface.  A  formulation  using  an  effective  transfer  speed 
associated  with  the  short  waves  gives  the  best  match  to 
the  data,  although  all  formulations  give  reasonable  re¬ 
sults.  Wave-resolved  conversion  of  the  wave-following 
measurements  to  a  fixed  reference  frame  show  that  the 
depth  penetration  of  the  turbulence  is  strongly  affected 
by  wave  orbital  motions,  such  that  turbulence  can  be 
carried  down  to  the  trough  level  and  below  (relative  to 
mean  sea  level).  The  depth  dependence  can  be  modeled 
as  a  power-law  decay  scaled  by  wave  height. 
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